Modern nuclear force predictions for the a-particle 
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We present new calculations of the a-particle which are based on the most modern nucleon-nucleon 
interactions alone and combined with the Tucson-Melbourne or the Urbana IX three-nucleon inter- 
action. Results for the binding energies and some properties of the wave function are given. On 
that phenomenological level little room is left for the action of a possible four-nucleon force. 

PACS numbers: 21.10.-h, 21.45. +v, 27.10. +h,21.30.-x 
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Few nucleon bound states have received increasing at- 
tention in recent years. The possibility of solving increas- 
ingly complex systems allows to probe the underlying dy- 
namics directly (for a detailed review, see fl). Several 
methods have been developed and applied to the 4N sys- 
tem and realistic forces, the GFMC §L the CHH j|, the 
SV @, the CRCGV § and the FY ||§ method. Very 
recently a new no core shell model calculation appeared 
H . The 4N system is an important test ground for both, 
the NN and the 3N nuclear interactions, because of its 
strong binding. In this article we want to address the 
question whether at least the ground state energy of the 
a-particle can be described by the most modern nuclear 
Hamiltonians. 

In recent years NN forces have been tuned very well to 
the rich NN data base, which led to a new generation of 
so called realistic NN forces: Nijm I, II [g], AV18 @ and 
CD-Bonn Q. While Nijm II and AV18 are purely local, 
Nijm I has a weak nonlocality in form of a V 2 -dependencc 
and CD-Bonn is quite nonlocal keeping the underlying 
Dirac structure of the one-boson-exchange without p/m 
expansions. All potentials describe the NN data base 
with a x 2 /data very close to 1. But they also have a 
large phenomenological character with a typical number 
of 40 fit parameters. The AV18 and CD-Bonn forces 
distinguish nn, pp and np interactions and thus include 
charge symmetry(CSB) and charge-independence(CIB) 
breaking. The interaction AV18 has in addition a whole 
set of electromagnetic corrections built in. 

It is well known that the two 3N bound states, 3 H and 
3 He, are theoretically underbound using only NN forces 
[jUI |l5|| . We compare our latest theoretical results based 
on fully converged Faddeev calculations to the experi- 
mental binding energies in Table |. The two nonlocal 
potentials lead to less underbinding than the local ones. 
The 3N calculations take the full CSB and CIB of the NN 
force into account including the isospin T — 3/2 admix- 
tures. In case of 3 He and 4 He the Coulomb interaction 
is included and, in addition, we use also the electromag- 
netic corrections in case of AV18. Because of an imple- 
mentation error the results for Nijm I, II changed with 
respect to |ll| . Table [| includes also the strongly model 
dependent kinetic energies which are correlated to the 



NN correlations. The local potentials have a somewhat 



harder core |13| leading to a higher kinetic energy. 

There are two additional dynamical ingredients, which 
should cure that underbinding, relativistic effects and 
three-nucleon forces (3NF). Though one can in principle 
always find an unitary (but very complex) transformation 
to a Hamiltonian without any 3NF Jl6| , we keep the given 
2N interactions for practical reasons and because of their 
generally accepted physical origin. In relation to those 
forces one or both of the two mentioned dynamical in- 
gredients are needed. There are still controversies about 
the role of relativistic effects [|l7]jl|] and as we think also 
open conceptual questions. Thus we consider here only 
a strictly nonrelativistic framework, having however in 
mind that relativistic corrections should finally be added. 

The topic of 3NF's is as old as nuclear physics |Tg| ] 
and based on meson exchanges various processes have 
been proposed in the past (for a review see pc| ). Among 
them the Fujita-Miyazawa force |2l| with an intermediate 
A generated by the exchange of two pions sticks out and 
is implemented in all modern 3NF models. Here we men- 
tion the rather popular 27r-exchange Tucson-Melbourne 
model (TM) p2j, the Brazilian version thereof |23| and 
the Urbana 3NF ||. There are also extensions to n- 
p and p-p exchanges [Q. The 27r-exchange model has 
been critically reviewed recently and a modified version, 
TM', has been proposed in |^5|,^6| which satisfies at least 
chiral symmetry. In this article we use the TM, TM' and 
the Urbana IX 3NF's. 

In the TM force enters the strong 7rNN vertex func- 
tion parameterized in the form of a monopole form fac- 
tor with a cut-off parameter A. Choosing it around a 
generally accepted value one achieves the right order of 
the lacking binding energy in the 3N system 1 27 14|,|2T 



Because there is a strong correlation between the a bind- 
ing energy, some 3N scattering observables |2!| and the 
3N binding energy, we fine tune A to the 3N binding en- 
ergies of 3 H and 3 He and this separately for each of the 
four NN forces. In this step we have not yet included 
T = 3/2 admixtures. We list the fit results in Table ||. 
Thereby the original TM parameters for the constants 
a,b,c and d have been used pfj| |. In case of the Nijmegen 
interactions we did not adjust A to the triton, because 
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E B [MeV] 



T [MeV] 



3 He 



E B [MeV] 



T [MeV] 



4 He 



E B [MeV] 



-8.012 
-7.623 
-7.736 
-7.654 



37.42 
46.73 
40.73 
47.51 



-7.272 
-6.924 
-7.085 
-7.012 



36.55 
45.68 
39.97 
46.62 



-26.26 
-24.28 
-24.98 
-24.56 



-8.48 



-7.72 



-28.30 



TABLE I. 3N and 4N binding energies for various NN potentials together with expectation values T of the kinetic energy. 



they do not include a modern specification of the 1 So nn 
force. Since for 3 H the lack of binding energy is a bit 
larger than for 3 He the A values for 3 H are slightly larger 
than for 3 He. 

Now having those 3N Hamiltonians at our disposal we 
can study the a-particle. For results referring to older 
forces see for instance |j],||j7j . We rewrite the Schrodinger 
equation into the Yakubovsky equations (YE) and 
thus decompose the wave function \& into 18 Yakubovsky 
components (YC). Due to the identity of the nucleons 
the YC's are not independent from each other and we 
can reduce their number to two: tpi an d ip2- Then the 
wave function reads 

* = (i - (i + p)p 34 )(i + P)ih + (l + P)(i + P)i>2 (i) 

and is expressed with the help of the permutations P — 
P12P23 + P13P23 and P = P13P24 where Pij are transpo- 
sitions of particles i and j. Thus P acts on the 3-body 
subcluster (123) and P interchanges the two two-body 
subclusters (12) and (34). This decomposition is highly 
advisable for scattering states since the boundary con- 
ditions can most easily be expressed for YC's. In our 
case of a bound state it would in principle be possible 
to solve directly the Schrodinger equation, but the usage 
of two YC's introduces in a natural manner two kinds of 
Jacobi coordinates which accelerates the convergence of 
a partial wave decomposition. 

For 4 identical particles the YE's reduce to two coupled 
integral equations 

i>i=G Q t 12 P [(I-P34) 1>i+W (2) 
V>2 = Go tu P [(1 - P34) i'i + H (3) 

Here in addition to the permutations the free 4N propa- 
gator Go and the NN t-operator ti2 occur. ti2 is driven 
by the NN force V12 through the Lippmann-Schwinger 
equation t 12 = Vi2 + V\i Go ti2- In case of 4N scattering 
one has to go one step further and solve the 3-body and 
2+2 subcluster problems beforehand in order to define 
the correct cut-structure ||. This is not necessary for 
bound states and the form of Eqs. (§)-@ is easier to 
handle numerically. One, presumably the most effective 
manner, to include 3NF's has been given in p^ ]. We stick 
to that. Then only the first of the two YE's is changed 
into 



*Ai = G t 12 P [(1 - P 3i )ipi + H + (1 + G fi2)GoV£gtf 

(4) 

(3) 

Here is that part of the 3NF which is symmetrical 
under exchange of particles 1 and 2. In case of the 27r- 
exchangc TM or the Urbana 3NF such a separation into 
3 parts is very natural. We solved the sets Eqs. (||)-(f|) or 
Eqs. (^)-@ in momentum space and in a partial wave 
representation. The first YC ip% stands for the "3+1" 
partition and is naturally described by two Jacobi mo- 
menta for the 3-body subcluster and one for the relative 
motion of the 4 th particle to the other 3. The second YC 
ip2 stands for the "2+2" partition and is naturally de- 
scribed by two relative momenta for the inner motion of 
the subclusters (12) and (34) and by the relative motion 
of the two subclusters. There are a lot of orbital and spin 
angular momenta as well as isospin quantum numbers to 
be coupled to J 71 " = + and T = 0. Both basis sets for the 
"3+1" and "2+2" partitions comprise about 1800 differ- 
ent combinations thereof in order to reach a converged 
description. The various momenta are discretized with 
roughly 35-45 grid points each. This leads to a huge ab- 
solutely full kernel matrix of dimensions 10 s x 10 8 . We 
solve the eigenvalue problem by a Lanczos type algorithm 
p3| , [34| and make intensive use of a massively parallel su- 
percomputer. 

We would like to remark that the introduction of the 
two kinds of Jacobi momenta in Eqs. (||),(||) and (^) 
leads to additional coordinate transformations which are 
hidden in the operator form of the YE's. They are equiv- 
alent to permutations. The crucial point in our calcula- 
tion is the treatment of these permutations and coordi- 
nate transformations. The direct interchange of arbitrary 
particles is unfeasible because of the huge dimension of 
the problem. It is therefore necessary to interchange par- 
ticles in two steps in such a way that at least one of the 
Jacobi momenta is not changed. This guarantees a block 
diagonal structure for the permutations. Only in this 
manner the calculation becomes feasible. For a detailed 
description see ||. 

Since we allow for CIB and CSB in the NN forces in 
principle the dominant total isospin state T — has to 
be supplemented by T — 1 and T = 2 admixtures. Our 
estimations lead to the result that their admixtures will 
change the binding energy only very slightly ( < 10 keV) 
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and thus at this stage we neglected them. But CIB and 
CSB lead to the prescriptions that the NN t-operators oc- 
cur in the form t = \ t n p + 3 t pp + i t nn in the NN isospin 
1 channels. This is different from the 3N system, where 
the corresponding linear combination is t = i t np + 1 t pp . 
The pp t-matrix also includes the effect of the Coulomb 
interaction. Since the bound nucleons are confined to a 
small space region, we can put to zero the Coulomb inter- 
action outside a radius of 10 — 20 fm. Then the Fourier 
transformation of this interaction is nonsingular. The 
results are cut-off independent and numerically stable. 

We show in Table | the a-particle binding energies us- 
ing NN forces only. As expected this theory underbinds 
the a-particle. The CD-Bonn result compares well with 
H . A correlation between the 3 H and 4 He binding en- 
ergies found in |36f| for simple forces remains valid also 
for the most modern ones. This is depicted in Fig. [j]. 
The experimental value is close to that straight line cor- 
relation, which nourishes the hope that by curing Bt one 
possibly will also cure B a . 
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FIG. 1. Correlation of 4 He against 3 H binding energies in 
MeV for the different potentials. The triangulars mark the 
predictions with 3NF from Table II. 

In Tabic |n] we present our results adding the TM 3NF 
adjusted individually to the 3N binding energies. We find 
a slight overbinding of about 300 keV for AV18, Nijm I 
and II . In case of CD Bonn the overbinding reaches 
about 800 keV. Note that adjusting the TM 3NF to 3 H 
leads to a somewhat larger overbinding. Nevertheless all 
these numbers indicate that with those Hamiltonians one 
reaches the a-particle binding energy rather closely and 
there is little room left for the action of 4N forces. On 
such a phenomenological level, however, it is not possible 
to decide about the need of 4N forces. It is always possi- 
ble to add another piece of a 3N force which in this case 
should be repulsive to reach the a-particle binding energy 
more accurately. For instance the ir-p exchange 3N force 
would provide at least one more parameter and both nu- 
clei, 3 He and 4 He, could be described. In that respect 
the Urbana 3NF is adjusted in such a manner jl| that 
in addition to the 3N binding energy also nuclear matter 
is taken into account. This then fixes a repulsive piece 



Model 


a-term 


6-term 


c-term 


d-term 


AV184 


-TM( a He) 


0.003 


-4.56 


-1.26 


-1.72 


AV184 


-TM'( 3 He) 


-0.26 


-3.90 


+3.00 


-1.02 



TABLE IV. Expectation values (in MeV) of the four parts 
of the TM 3NF model with respect to wave functions gener- 
ated with AV18+TM and AV18+TM' for 4 He. 

in that 3NF, with the result that the a-particle binding 
comes out essentially right. We repeated that calculation 
for the Urbana 3NF firstly performed with the GFMC 
method gj. Our result based on YE's is given in Table [] 
together with the previous one. There we also show a 
triton result by the CHH-method p5[ . We estimate our 
numerical accuracy to be about ±3keV (±50keV) for the 
3N (4N) system. There appear to be small differences 
between our and the GFMC results, especially in the ki- 
netic energies. 

Like in 3 H one can also separate the 4 He wave func- 
tion into S, S', P and D-state probabilities. Here S is 
spatially symmetric, S' has two-dimensional mixed sym- 
metry and P and D are the total L = 1 and 2 orbital 
angular momentum parts . We compare the two nuclei in 
Tables ILL For 4 He S' is reduced, P somewhat enhanced 
and probabilities shifted from S to D. It is remarkable 



that the TM 3NF together with CD-Bonn reduces the 
D-state probability, in contrast to all other cases. 

In it is argued that chiral symmetry requires that 
the c-term in the TM 3NF should be dropped, leading 
to a TM' 3NF (one keeps the remaining constants b,d 
unchanged and a is replace by a' = a — 2c). As seen 
in Table [b] the resulting a-article binding energy after 
fitting the cut-off to 3 He coincides essentially with the a- 
particle binding energy. It is also interesting to see that 
the c-term has a significant effect on the wave function 
as demonstrated in Table |y|. There we show the expec- 
tation values of the four different terms contributing to 
the TM 3NF evaluated with the wave function including 
TM and the a, b and d expectation values based on the 
wave function including TM'. In relation to that latter 
wave function one can also evaluate the c-term. Inter- 
estingly it turns out to be repulsive, whereas for the TM 
wave function it is attractive. Also it is interesting to see 
that for TM the a-term is negligible, whereas it has some 
importance for TM'. These examples demonstrate again 
that 3NF's change wave functions and cannot be treated 
perturbatively, a fact known since long time ]2^] . 

Summarizing, after adjusting 3NF's to 3N binding en- 
ergies the a-particle binding energies based on modern 
nuclear forces are rather close to the experimental value. 
This indicates that 4N forces (at least for T = 0) will 
be unimportant. More details on wave function proper- 
ties will be published elsewhere. The FY equations are 
perfectly well under control for realistic NN and 3NF's 
and will be a perfect tool to study upcoming new force 



3 



Potential 


A [m*] 


3 H 




He 


4 


He 






E B [MeV] 


T [MeV] 


E B [MeV] 


T [MeV] 


E B [MeV] 


T [MeV] 


CD Bonn + TM 


4.784 


-8.480 


39.10 


-7.734 


38.24 


-29.15 


83.92 


CD Bonn + TM 


4.767 


-8.464 


39.03 


-7.720 


38.18 


-29.06 


83.71 


AV18 + TM 


5.156 


-8.476 


50.76 


-7.756 


49.69 


-28.84 


111.84 


AV18 +TM 


5.109 


-8.426 


50.51 


-7.709 


49.47 


-28.56 


110.92 


AV18 + TM' 


4.756 


-8.444 


50.55 


-7.728 


49.54 


-28.36 


110.14 


Nijm I + TM 


5.035 


-8.392 


43.35 


-7.720 


42.59 


-28.60 


93.58 


Nijm II + TM 


4.975 


-8.386 


51.02 


-7.720 


50.13 


-28.54 


113.09 


AV18 + Urb IX 




-8.478 


51.28 


-7.760 


50.23 


-28.50 


113.21 


AV18+Urb IX (GFMC)J| 




-8.47(1) 


50.0(8) 






-28.30(2) 


112.1(8) 


AV18+Urb IX (CHH) pf 




-8.476 


51.26 










Exp. 




-8.48 




-7.72 




-28.30 





TABLE II. Cut-off parameters A, adjusted 3N binding energies, resulting a particle binding energies for various force 
combinations. Expectation values of the kinetic energy are also shown. Bold faced results have been adjusted to the experiment. 



Model 


S [%] 


S' [%] 


4 He 

P[%] 


D [%] 


S [%] 


S' [%] 


3 He 

P[%] 


D [%] 


AV18 


85.45 


0.44 


0.36 


13.74 


89.95 


1.52 


0.06 


8.46 


AV18+TM( 3 He) 


85.10 


0.30 


0.75 


13.84 


89.86 


1.26 


0.15 


8.72 


AV18+TM'( 3 He) 


83.27 


0.31 


0.75 


15.68 


89.46 


1.25 


0.13 


9.16 


AV18+Urb-IX 


82.93 


0.28 


0.75 


16.04 


89.39 


1.23 


0.13 


9.25 


CD-Bonn 


88.54 


0.50 


0.23 


10.73 


91.45 


1.53 


0.05 


6.98 


CD-Bonn+TM( 3 He) 


89.23 


0.43 


0.45 


9.89 


91.57 


1.40 


0.10 


6.93 



TABLE III. S, S', P and D state probabilities for a and 3 He 



structures given in chiral perturbation theory |37]|38|] . 
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